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Abstract 

We discuss the analytic properties of the Callan-Symanzik /3-function ()[g) associated 
with the zero-momentum four-point coupling g in the two-dimensional (f/^ model with 
0{N) symmetry. Using renormalization-group arguments, we derive the asymptotic be- 
havior of I3{g) at the fixed point g* . We argue that (3'{g) = (3'{g*) + 0{\g - g*\^^'^) for 
iV = 1 and P'{g) = (3'{g*) + 0(1/ log \g - g*\) for > 3. Our claim is supported by 
an explicit calculation in the Ising lattice model and by a 1/A^ calculation for the two- 
dimensional (j)^ theory. We discuss how these nonanalytic corrections may give rise to a 
slow convergence of the perturbative expansion in powers of g. 



1 Introduction 



Renormalization-group theory is a very important tool for the understanding of the critical be- 
havior of statistical models in the neighbourhood of the critical point. We consider models with 
an A^- vector real order parameter and 0{N) symmetry. Because of universality, quantitative 
predictions can be obtained by studying any theory belonging to the same universality class. 
For the models we are dealing with here, we may consider the Ginzburg-Landau Hamiltonian 
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where is an A^-component real field. This Hamiltonian describes many interesting systems 
at criticality. The liquid-vapour transition in fluids and the infinite-length properties of poly- 
mers in dilute solutions correspond to the N = 1 (Ising) and = model respectively; the 
^He superfiuid phase transition is in the same universality class of the three-dimensional two- 
component theory {XY model), while the Hamiltonian (|l]) with = 3 describes isotropic 
ferromagnetic materials. Three-dimensional A^-vector systems and two-dimensional systems 
with N < 2 have a conventional critical behavior: Thermodynamic quantities have power-law 
singularities near the critical point. On the other hand, in two dimensions the XY model shows 
a Kosterlitz-Thouless transition, while for A^ > 3 no finite-temperature transition exists: The 
correlation length diverges only for T — > 0. For A^ > 3 the theory is asymptotically free with a 
critical behaviour described by the perturbative renormalization group applied to the nonlinear 
cr-model. 

Precise estimates of the critical parameters in the symmetric phase can be obtained using 
several different methods. One of them, which provides in many cases very precise results, 
relies on a perturbative expansion in powers of the zero-momentum four-point renormalized 
coupling g performed at fixed dimension d The theory is renormalized by introducing a 
set of zero-momentum conditions for the (one-particle irreducible) two-point and four-point 
correlation functions: 

r^'\p)ap = 5ap [m'+p' + 0{p')] , (2) 

r*^^) (0, 0, 0, 0)af3'r5 = Z(.^ m^""^ g ^ (dafsS^S + Sa-ySfSS + SasSfB-,) ■ (3) 

For m — > 0, the coupling g is driven toward an infrared-stable zero g* of the corresponding 
Callan-Symanzik /3-function 

(4) 

The derivative of the /5-function at g*, (3'{g*), is related to the leading nonanalytic correction-to- 
scaling exponent. Usually — but we shall argue here that this may not always be the case — the 
leading nonanalytic corrections are determined by the critical dimension uJi of the leading 
irrelevant operator: in this case, we have 13' {g*) = uji. At present, (3{g) has been computed to 
six loops in three dimensions and to five loops in two dimensions 

Perturbative expansions in powers of g are asymptotic. In order to obtain estimates of 
universal critical quantities, it is essential to resum the perturbative series. This can be done 
by exploiting their Borel summability and the knowledge of their large-order behavior (see 
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e.g. and references therein). The large-order behavior of the series S{g) = J^^kO^ is related 
to the singularity gi, of the Borel transform B{g) that is closest to the origin. For large k, 



with 



(5) 



The value of gh can be obtained by means of a steepest-descent calculation P,^]. It depends 
only on the Hamiltonian, while the exponent h depends on which Green's function is considered. 
If the perturbative expansion is Borel summable, then gi, is negative. Since the Borel transform 
is singular for g = gt, its expansion in powers of g converges only for \g\ < \gb\. An analytic 
extension can be obtained by a conformal mapping [0, such as 



y(9) 
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1 - g/gb + 1 



(6) 



The Borel transform becomes an expansion in powers of y{g) that converges for all positive 
values of g, provided that all singularities of the Borel transform are on the real negative 
axis . Therefore, the use of the Borel transform and of the conformal mapping (^) transforms 
the original asymptotic series into a convergent expansion. Any universal quantity, such as 
the critical exponents, is estimated by resumming the corresponding perturbative series and by 
evaluating the resummed function of g at the fixed-point value g* . 

The critical value g* of the renormalized coupling is a universal quantity. Therefore, it can 
also be obtained by considering any statistical (lattice) model belonging to the corresponding 
universality class. Then 

' SAT X4 " 
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where t = T/Tc — 1, x is the magnetic susceptibility, ^ the second-moment correlation length, 
and X4 the zero-momentum four-point connected correlation function. Using Eq. (|^, one can 
obtain an independent estimate of g*. 

An important issue in the field-theoretical (FT) approach concerns the analytic properties 
of P{g). General renormalization-group arguments (see also |TD|,|TT[) and explicit calcula- 

tions to next-to-leading order in the framework of the 1/A^ expansion [|1^,|13[ show that (3{g) is 
not analytic aX g = g* . This fact may cause a slow convergence of the resummations of the per- 
turbative series to the correct fixed-point value. The reason is that this resummation method 
approximates the /3-function in the interval [0, g*] with a sum of analytic functions. Since, for 
g = g*, the /3-function is not analytic, the convergence at the endpoint of the interval is slow. 
This may also lead to an underestimate of the uncertainty that is usually derived from stability 
criteria. In spite of these problems, in three dimensions, FT results are in good agreement]^ 

^ Small discrepancies are only observed for = and = 1. For instance, we may compare the estimates 
of g* and lui obtained using the fixed-dimension FT approach with the apparently best estimates obtained from 
the analysis of high-temperature (HT) expansions and from Monte Carlo (MC) simulations for lattice models 
in the same universality class. For = 1, the analysis of the fixed-dimension FT expansion gives g* = 1.411(4) 
and ui = 0.799(11) (m), to be compared with the lattice resuhs g* = 1.402(2) ||l|l (HT) and ui = 0.845(10) |l6| 
(MC). The results are in better agreement for N = 2: the analysis of the fixed-dimension ^-expansion leads to 
g* = 1.403(3) and wi = 0.789(11) (l|, to be compared with g* = 1.396(4) (HT) and wi = 0.79(2) |l| 
(MC). 
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with the estimates obtained in other approaches 1 12, 14^19], showing that the above-mentioned 
nonanalyticity causes only very small effects that are negligible in most cases. Using general 
renormalization-group arguments, for three-dimensional models one expects |Q| 



m = -P\9*W ~9)[l + - gf + a^{9* - g) + ■ ■ ■] . (8) 

where (3'{g*) = ui and p is a noninteger exponent that is equal to the smallest of the following 
exponent combinations: p = UJ2/UJ1 — 1 where UJ2 is the scaling dimension of the next-to-leading 
irrelevant operator, p = 1/A where A = ujiu, and p = 7/A — 1. Notice the last exponent that 
was neglected in P,|12| and that is due to a subleading correction in g(t) proportional to f^. 
Such a term is related to the presence of an analytic background in the free energy. For small 
values of A^, we have A = ujiu ^ 1/2, 002/001 ~ 2 pO|, and 7/A > 2, so that p = UJ2/0J1 — 1 ~ 1. 
In this case the leading nonanalytic term is practically undistinguishable from the analytic 
one, and therefore, one expects only small systematic deviations. For increasing values of A^, 
p decreases, but at the same time Oi — > 0. Thus, also in this case we expect the nonanalytic 
terms to give rise to small systematic deviations. 

The situation worsens in the two-dimensional case which we consider in this report. As a 
matter of fact, at variance with the three-dimensional case, two-dimensional FT estimates are 
much more imprecise [Q]. We shall argue here that the large observed deviations are caused 
by the nonanalyticity of the renormalization-group functions at g* . In order to support this 
argument, we shall compute the behaviour of the /3-function for g g* m two cases in which 
exact results can be obtained exploiting different techniques. 

First, we shall address the A^ = 1 case (i.e. the Ising universality class) in which conformal 
field theory (CFT) techniques allow the determination of the whole spectrum of relevant and 
irrelevant operators of the theory. We shall first show that CFT predicts tui = 2 for the 
renormalization-group dimension of the leading irrelevant operator, excluding oji = 4/3, as 
it has been claimed sometimes. Then, we will consider the lattice Ising model and we will 
show that Eq. (P) holds with P'{g*) = 7/1/ = 7/4 and p = 1/7. Notice that in this case 
P'{g*) 7^ t^i = 2. We will then argue that this is the generic behaviour one should expect for 
models in the Ising universality class. At variance with the three-dimensional case, here p is 
very small and thus it may be responsible for large systematic deviations in the resummation 
of the perturbative series. In App. B we study some simple Borel-summable asymptotic series 
behaving as (H) with p = 1/7. We apply the resummation method describe above, finding very 
poor estimates of (3'{g*) with largely underestimated error bars. 

Second, we shall study the multicomponent (f)^ theory with A^ > 3. Since the model is 
asymptotically free, we can predict oJi = 2 and we can show that logarithmic corrections should 
be expected at the critical point. A large- A^ calculation confirms the theoretical predictions. 

2 N =1 (j)^ theory in = 2 

Let us first consider the Ising case, i.e. the case in which the field in the 0^ Hamiltonian is 
a one-component real field. In the four-loop series of (3{g) is analyzed using the resummation 

^ This is due to the fact that the correction-to-scahng term with the smallest exponent appearing in g{t) is 
V and not f^^^ . A correction term proportional to in g{t) is due to the presence of an analytic background 
in the free energy. 
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procedure presented in the introduction: they obtain g* = 15.5(8) and P'{g*) = 1.3(2). Refer- 
ence PI computes the five-loop contribution and presents an analysis of the extended series using 
a Pade-Borel resummation: they obtain^ g* = 15.39(25) and P'{g*) = 1.31(3). These results 
for g* do not agree with the very precise estimates obtained by a transfer-matrix analysis of the 
standard square-lattice Ising model g* = 14.69735(3), and by exploiting the form-factor 
bootstrap approach |Q, g* = 14.6975(1) (see also ||12|,^,|2^ for high-temperature results). 



The result for (3'{g*) has been interpreted |]^, ^ as an indication in favor of the exact result 
(3'{g*) = 4/3 that would imply the existence of an irrelevant operator with uJi = 4/3. However, 
the corresponding scaling corrections do not appear in the standard lattice Ising model in which, 
thanks to the known exact results (see e.g. p5|^7|), a detailed analysis of the leading correction 



terms is possible. In principle, this fact does not imply that the interpretation of is wrong, 
since it could be simply explained by the absence of the corresponding irrelevant operator in 
the lattice Ising model, which is only one of the possible realizations of the 0^ universality class. 
However, we shall show below that this is not the case and that no subleading operator with 
uji =4/3 exists in any unitary model belonging to the Ising universality class. In particular, it 
does not exist in the N = 1 (p^ theory. 

Let us briefiy comment on this last point. The ui =4/3 interpretation was supported by 
the fact that an operator with renormalization-group dimension ui = 4/3 exists in a particular 
nonunitary extension of the Ising universality class which is conjectured to describe Ising perco- 
lation. However, such an operator can only exist in nonunitary theories, and as a consequence, 
it cannot be observed in the unitary 0^ theory. We shall argue in this paper that the estimate 
of u!i obtained in the framework of the perturbative expansion at fixed dimension is strongly 
affected by nonanalytic corrections in the /3-function. The fact that one obtains (3'{g*) ~ 4/3 is 
only a coincidence and is not related to the presence of the nonunitary operator with ui = 4/3 
mentioned above. In order to clarify the issue we have added in App. ^ a discussion on the 
nonunitary extension of the Ising universality class and its relation with the Ising percolation 
problem. 

The only ingredients that are needed to extend the Ising result — the absence of an exponent 
Ui = 4/3 — to the most general unitary model in the N = 1 (p^ universality class are Wilson's 
renormalization group and some basic results of CFT. 

In Wilson's approach, we can rewrite 7i as 

n = n* + Y.uo{m)o, (9) 

{O} 

where H* is the fixed-point Hamiltonian, {O} a complete set of operators, and uoijn) the 
corresponding nonlinear scaling fields depending on the inverse correlation length m. Then, we 
observe that the 0^ theory is unitary. This can be proved to all orders of perturbation theory. 
It can also be proved nonperturbatively by considering the lattice regularization of the model 
d^). Indeed, the lattice theory corresponding to (|I|) — and, of course, also the standard Ising 
model which is a particular limit of the lattice 0^ theory — with nearest-neighbor couplings is 
exactly refiection positive, a property that guarantees the unitarity of the Minkowski theory. 
At the critical point the theory becomes conformally invariant. Now the main point is that 



•^We applied the Le Guillou-Zinn- Justin resummation method 1^,0], using the conformal mapping (^), to the 
five-loop series of [p|: we obtained substantially equivalent results. 
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in the framework of CFT there exists a complete classification of all possible Z2 symmetric 
unitary theories [^,|2^. Moreover, their operator content is exactly known. This means that 



all dimensions of the operators O that may appear in Eq. giving rise to a unitary theory 
are exactly known. ^ In particular, no operator with dimension cui = 4/3 exists. 



According to the CFT analysis ||21], the leading irrelevant operator is TT, where T 
denotes the energy- momentum tensor, which is expected to give rise to corrections of order t^, 
t being the reduced temperature. On the square lattice — but not in a rotationally-invariant 
model or on lattices with different rotational symmetry, for instance, on the triangular lattice — 
one must also consider a second operator, T = + T^, which is degenerate with the first 
one. While TT is rotationally invariant, T breaks rotational invariance and has only the 
reduced symmetry of the square lattice. Since correlation function of T with rotationally 
invariant operators vanish, such operator should not contribute at order to observables that 
are rotationally invariant, but only at order (indeed, (7^7^(9) does not vanish even if O is 
rotationally invariant). Of course, T should contribute to order to observables that have an 
angular dependence (an explicit example will be given below). 

In the last years there has been extensive work trying to understand the origin of the 
subleading corrections in the lattice Ising model. The unexpected result is the fact that no 
correction-to-scaling term due to TT has been observed. Let us review the evidence for this 
fact: 



1. The analysis of the susceptibility [p2H34|l for /i = indicates that the corrections of order 
t, t^, can be interpreted as purely analytic ones. 



2. The analysis of the free energy on the critical isotherm as a function of h [31 1 does not 



find any evidence of correction-to-scaling terms that can be associated to TT . 

The analysis of the free energy, correlation length and susceptibility at the critical point 



in a finite box [^, Q shows the presence of corrections with ujx = 2. These corrections 
however appear to be due to T only. Indeed, they are not present on the triangular and 
honeycomb lattices — on these lattices T cannot contribute and the first expected 
correction has ui = 4 — and moreover, the dependence of these corrections on the shape 
of the box is consistent with the behaviour expected for a spin- four operator as T is ||36||. 



Here, we want to add further evidence for the absence of TT by considering the observables 
characterizing the large-distance behaviour of the two-point function on a square lattice. Indeed, 



for a; —> 00 we can write [37 



where = 4 X^/^ sin^(p^/2) and the integration is extended over the first Brillouin zone. The 
quantities Z{P) and M{P) are known exactly [^. For t = 1 — /?//5c — 0, we can write 



Z{(3) = (l28V2py^\y%l\l + 0{ut)], (11) 



"^Let us stress that our argument is by no means original. It was, for instance, already present in that 
appeared right after the classification of unitary CFT's. In this respect, our main new contribution is the exact 
calculation of p and the use of this result (discussed in detail in App. ^) to show the relevance of nonanalytic 
corrections in the FT Callan-Symanzik /3-function. 
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(12) 



where Ut is the nonhnear scahng field associated with the reduced temperature at zero magnetic 
field h and Vh is related to the nonlinear scaling field Uh associated with hhj Uh = hv^ + 0{h^). 
Exphcitly [|3,|3|,||] 



Vh 



\ V2 6 



1 + + 

72 16 



6^2 



(13) 
(14) 



Using (p!oD we can derive the angle-dependent correlation length ^{6) defined from the large- 
distance behavior of the two-point function along a direction forming an angle 9 with the side of 

we obtain! 



the lattice. Using the expression of ^{6) in terms of M(/5) reported, e.g., in [0,^ 



4:PcUt 



l + i/3^os(40)«? + O« 



(15) 



Thus, we see analytically that no correction of order 0{t^) appears in the on-shell renormal- 
ization constant Z{(3) — both TT and T are absent. In ^{6) a 0{t'^) correction does appear as 
already observed in [^. However, it is proportional to cos(4^), and thus it is due only to the 
leading operator breaking rotational invariance. No contribution from the rotationally-invariant 
operator TT appears. 

The expansion of the Callan-Symanzik /5-function can be derived using the same arguments 
employed by Nickel in three dimensions. Let us first consider the lattice Ising model and 

the coupling g{t) defined in (|^ as a function of the reduced temperature. The expansion of x 
and X4 is well established: 



X 

X4 



C2Ut 



7/4 2 
Vh 



7/4 
PlUt 



11/4, , 11/4 , 

P2Ut log Ut + PsUt + 



^ -xi/^ 4 /-I , 11/4 , 



(16) 
(17) 



where C2, C4, pi, p2, Ps, and are known constants [^T], |2^, ^ |33|, ^ . In particular 
Pi = —0.1081812 . . . Next, we determine the asymptotic behavior of fj,2 = J2x ^'^{o'o'^x) = 4X'C^ 
from its high-temperature expansion (HT). The analysis of the 52nd-order HT expansion^ of 
/i2 shows that its Wegner expansion can be written as 

/i2 = A2U;'''^\l (1 + + ...). (18) 

The constant ps has been computed with high accuracy in the following way. We have first 
defined a new series s obtained by expanding in powers of /? the quantity {n2ul^ v'^'^ / A2 — l)u'i~'^ , 
where A2 = 1.238136098, and Ut, Vh are given by Eqs. (|T3|) and (|1^) truncated at order 



^In particular, the correlation lengths along the side {6 ~ 0) and the diagonal {9 — it /A) of the lattice are 
respectively given by = — lntanh/3 — 2(3 and = — \/21nsinh2/?. 

^ The HT expansion of ^2 can be found to 0{I3'^^) in Ref. |^. The 52nd-order series has been kindly provided 
by Tony Guttmann 
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included. Then, we analyzed s by means of first-order inhomogeneous integral approximants 
biased to have a singularity at /3 = f3c- We verified that the critical exponent associated to 
the singularity is positive and then computed the value of s for /3 = Pc- We obtain finally[] 
P5 = -0.388720(3). It follows that 



9{t) = g*\l- PiuJ^ - P5uf + 0{u^'^ logut) 



Since the second-moment mass m(t) = l/^(t) = (4x//U2)^''^ scales as 



m(t) 



2 _ 4<^2 2 



1 + 0(u 



7/4^ 



(19) 



(20) 



we obtain for the square-lattice Isin^ /3-function 



i3{g) = m 



dg 



2m^ 



'dm? 



dg 
dut 
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/\g[l + h\^g\'l' + h2\Ag\"' + fesi A^7|^/^ + ■ 



2/7 



1 3/7 



c?m \ (iMf ) dut 4 

(21) 

where A(yf = g* — g, and for the nonuniversal constant 6i, 6i = p^{—g*Pi) ^^''/(7pi) ~ 0.480(4). 
It follows that P'{g*) = 7/4 and p = 1/7. Let us stress again that this value of I3'{g*) is 
not related to the exponent of the leading irrelevant operator that we expect to be two. This 
phenomenon occurs whenever 7 < ujiv. Indeed, in g{t) there is a correction-to-scaling term 
proportional to because of the presence of an analytic background in the free energy ||32|| . 
If 7 < ciJiZ/, it represents the leading nonanalytic correction in g{t) and therefore (3'{g*) = 
^jv ^ It should be noted that such a phenomenon does not arise in three-dimensional 
0{N) models, where the leading analytic corrections are determined by the leading irrelevant 
operator. For instance, for the three-dimensional Ising model 7/1/ = 2 — ^ 1.96 > cji ~ 0.8. 
We also mention the recent result (3'{g*) — 1.88 obtained in HBI using a numerical approach 
based on the high-temperature expansion of the Ising model, which is not too far from our 
exact prediction 7/4. 

Now, the question is: which behavior should we expect for the 0^ field theory? In other 
words, does Eq. ([19D holds for a generic model in the N = 1 (p^ universality class or are some 
terms absent? And, in particular, are the conditions Pi ^ and P2 7^ a particular feature of 
the lattice Ising model only? Sometimes, see e.g. |^ and the discussion of ||12|, it is conjectured 



that the /3-function is analytic in FT models. However, it was shown in [12| that this conjecture 
is not true: In the large-A^ limit, nonanalytic terms are indeed present. Unfortunately, in the 



^ It is worth noting that ps w — = —0.388722... within error bars, so that fi2 can be written as 
^2 = A2vf^Uf ^^^^ — + • ■ • This equation may be explained in terms of a momentum-dependence of the scahng 
field Uh- Indeed, fj,2 is not a zero-momentum quantity and thus it is related to the free energy in the presence 
of a nonuniform magnetic field h(x). But, in this case we expect additional contributions to the scaling fields, 
proportional to derivatives of h{x) [Q. Our result for /i2 can be explained if the scaling field Uh is a functional 
of h{x) with a small-momentum behavior Uh — Uh\h=const — ^d'^h{x)+ higher derivatives. Note also that the 
term in ^2 cannot be interpreted as a contribution due to irrelevant operators. Indeed we do not expect 0{u1) 
contributions associated with TT, nor with the nonrotationally invariant T, since ^2 is a rotationally invariant 
quantity. This point needs further investigation. 



^ Notice that in the lattice model g approaches g* from above as t 
happens. For a discussion see |45| and references therein. 



0, while in the FT model the opposite 
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two-dimensional case for = 1, we do not have any analytic control on the corrections to 
l3{g). Nonetheless, we conjecture that Eq. (0) holds also for the FT N = 1 model — of course, 
with different coefficients bi, 62 since the /3-function is not universal. We have essentially two 
arguments to support our conjecture: 



a] We do not see any reason why the bulk term that originates the pit'/'^ contribution in ([19D 

should not be present. Indeed, the analytic contribution is not a lattice artifact but has a 
well-defined FT meaning. In the CFT framework, it can be considered as a signature of 
the Identity operator and of its conformal family. Thus, also for the FT model, we expect 

Pi 7^0. 

b] A correction is certainly present in g{t), since we expect the operator TT to be present 

in the FT model. Thus will not be zero in (|1^), although it will be no longer related 
to the correction appearing in ^2- 

It is important to note that the strong nonanalytic corrections at g = g* we have found may 
explain the large observed deviations among the perturbative FT estimates of g* and (5'{g*), 
the high-precision numerical results for g*, and our prediction for (3'{g*). As a test, in App. 
we have considered a simple Borel-summable function that has an asymptotic behavior of the 



form (pTI). We have applied the standard resummation method presented above, observing large 
systematic deviations at g = g* and a systematic underestimate of the error bars. We should 
note that these discrepances, although provide support for the presence of strong nonanalytic 
corrections at g = g*, do not support our specific expansion (|2l|) . Indeed, even if pi = in 
(|T9|), neglecting logarithmic terms, we would obtain 



Pig) = -2Ag[l+c,\Ag\'/' + ...). (22) 
Thus, also in this case, there would be a strong nonanalytic correction. 



3 N >3 (j)^ theory in d = 2 



Let us now consider the multicomponent 0^ theory with A^ > 3. For A^ = 3, the Pade-Borel 
analysis of the five-loop series yields the estimates g* = 12.00(14) and P'{g*) = 1.33(2). The 
result for g* is in reasonable agreement with the more precise estimate g* = 12.19(3) obtained 
by employing the form-factor bootstrap approach [^,|^. We shall now argue that the estimate 
j3'{g*) ~ 4/3 is again incorrect and that the correct value should instead be P\g*) = 2. 

The standard scenario predicts that, for A^ > 3, the theory is massive for all temperatures. 
The critical behavior is controlled by the zero-temperature Gaussian point and can be stud- 
ied in perturbation theory in the corresponding A^-vector model. One finds only logarithmic 
corrections to the purely Gaussian behavior. It follows that the operators have dimensions 
that coincide with their naive (engineering) dimensions, apart from logarithmic multiplicative 
corrections related to the so-called anomalous dimensions. The leading irrelevant operator has 
dimension two and thus, for m — 0, we expect |^ 



g{m) 



9* {I 



cm 



^ ( — In m 



2^< 



1 + 



'ln(— Inm^ 
Inm^ 



(23) 
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where C is an exponent related to the anomalous dimension of the leading irrelevant operator, 
and c is a constant. A one- loop calculation in the framework of the O(A^) cr model gives 
( = 2/{N — 2) [^. Differentiating with respect to the mass, one obtains 



^to)^™^^-2As(l + i^ + ...). (24) 

with Ag = g* — g. Therefore, one expects I3'{g*) = 2 with logarithmic corrections. 



The expansion ( p^) for (3(g) is confirmed by a next-to-leading order calculation in the 
framework of the large- expansion. Indeed, using the expression for i3{g) reported in |]l3l and 
performing an asymptotic expansion around g* (see App. C for details), one finds 



Pis) = -2(9- 




Ine I Iney 21n20 Un^Q 



(25) 



where /(9) = ln(— 21n9) and Q = {g* — g)/g*. Comparing Eq. (^ with Eq. (0) we obtain 
( = 2/N + 0{1/N'^), in agreement with the above-mentioned result ( = 2/{N — 2). 

Thus, for A^ > 3 we predict very strong nonanalytic corrections aX g = g*. A numerical 
study on a function with the asymptotic behavior (p^) (see App. B) shows that such corrections 
give rise to a slow convergence of the perturbative resummations. In particular, the estimate of 
Ui may be incorrect in spite of the stability of the results with the number of loops considered 
in the analysis. It is thus not surprising that Ref. |^ find l3'{g*) ~ 4/3 instead of the correct 
result I3'{g*) = 2. 
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A Nonunitary extension of the Ising model 

The 4/3 operator appears in a nonunitary extension of the Ising model that describes Ising 
percolation. 

Let us first of all explain what we mean with the notion of "nonunitary extension" of the 
Ising universality class. The starting point is the classification of the minimal unitary conformal 
field theories discussed in [p^,|29[ . 

The operator content of the unitary CFT's that only possess a symmetry (like the 0^ 
theory and its multicritical generalizations) is defined by the weights: 

^'^ 4m(m + 1) ^ ^ 
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with m = 3, 4, 5 • ■ ■ and the constraints 1 < p < (m — 1), 1 < q < p. The relation between h 
and the renormahzation-group eigenvalue y is y = 2 — 2h. For the Ising model m = 3. Higher 
values of m correspond to multicritical Ising-like models (i.e. theories with a Z2 symmetric 
potential with powers up to These are the continuum-limit CFT's that correspond to 

the models introduced in |^T],0. With m = 3 we have only three allowed combinations of 
{p,q): (1,1)) (2,1) and (2,2) that correspond to the identity, energy and spin operators of the 
Ising model. They are called "primary" fields. From any one of these primary fields one has 
then an infinite tower of "secondary" fields whose scaling dimensions are shifted by integers 
with respect to those of the primary fields. Since in the Ising model all the primary fields are 
relevant, all the irrelevant fields must be shifted by integers, hence they cannot be distinguished 
from the analytic corrections. This is the only model in which this happens. In all other models 
with m > 3, there are primary fields that are irrelevant and hence are candidates for nontrivial 
subleading scaling dimensions. 

Besides unitary theories, there is an infinite set of nonunitary ones for all the rational (but 
noninteger) values of m. Apart from the fact that they do not fulfill unitarity, they have the 
same properties of those with integer m. In particular, their operator content is completely 
known and closed expressions for the correlators exist. These models (with both integer and 
non integer values of m) are usually called Rational Conformal Field Theories (RCFT). 

However this is not the end of the story. In the last few years it has been realized that it 
is possible to give a meaning, in the framework of the so called Logarithmic Conformal Field 
Theories (LCFT) |]53|, also to more general theories, obtained by including in the operator 



algebra some of the operators corresponding to the values of p and q excluded in eq. ( [A.l|) []54 



For instance, in the Ising case (i.e. m = 3) in which we are interested one should enlarge 
the set of operators of the standard Ising CFT to those of the type /i3_„, n = 1, 2, ■ ■ ■ and /ifc^4, 
with k = 1,2 ■ ■ ■. The LCFT obtained in this way is what we mean by "nonunitary extension" 
of the Ising model. 

Despite the fact that these theories are much more difficult to study than the standard 
RCFT's, several interesting results have been obtained in these last years (for a recent account 
see for instance ||5^,^ and references therein). For the purpose of the present paper we only 
need to know the scaling dimensions of the new operators. These can be easily obtained by 
looking at eq. (|A.1|) . 

In particular, in the Ising case, we see that h^^i = | hence 1/3,1 = — |, which is exactly 
the irrelevant operator that we are looking for. Further examples of such operators (only the 
relevant ones are listed) are: 

h3,2 = § hence 2/3,2 = ^, 

/i3,3 = I hence ?/3,3 = |, 

^2,4 = jQ hence ?/2,4 = ^• 

Note that, for all the values m > 3 (i.e in the multicritical models), Eq. (|A.1| ) admits a 
unitary, well defined, operator of type h^^i with weight (m + 2)/m so that y = — 4/m. Thus, a 
naive limit m ^ 3 would lead to an operator with y = —4/3, This argument is usually given 
to support the existence of a scaling operator with cui = 4/3 (see e.g. 0]). However, as we have 
seen, exactly for m = 3 this operator becomes "border-line" and it does not belong anymore 
to the Ising universality class, but only to its nonunitary extension. Thus, the limit m — > 3 of 
Eq. ( [A.l|) cannot be considered as an indication in favor of the presence of a tui = 4/3 field in 
the (unitary) Ising universality class, which the (f)^ theory belongs to. 
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Another context in which the y = 4/3 field appears, which is completely independent and 
allows to share some more light on its meaning, is the Coulomb gas approach to the g-state 
Potts models due to Nienhuis |5^. By mapping the Potts model in a suitable vertex- type 



model Nienhuis was able to identify both the leading and the subleading thermal and magnetic 
operators as a function of q. For q = 2 the subleading thermal operator is exactly y = —4/3 



and the subleading magnetic operator is y = 13/24 (see also [58|). However, as already noted 
in |B^, these are operators of the Vertex model and not of the Ising model and they decouple for 



q = 2. In other words, the Vertex model of Nienhuis is a good candidate for an exactly solvable 
model whose continuum limit is the nonunitary extension of the Ising model. If one requires the 
Vertex model to have a "physical spectrum" according to the definition given in |^ , then one 
selects only the operators of the standard Ising model and the y = 4/3 operator decouples. The 
requirement of having a "physical spectrum" is equivalent to impose unitarity on the model. 

It would be nice to have some kind of insight of the physical meaning of the above-mentioned 
operators directly from Ising model. Some hints in this direction are given by the so called "Ising 
percolation" problem i.e. the behaviour of the Coniglio-Klein clusters in the Ising model. It 
turns out that the relevant operators in the nonunitary extension of the Ising universality class 
(i.e. both the standard ones y = 1 and y = 15/8 and the "border-line" ones y = 13/24, 
y = 11/8 and y = 5/3) become fractal dimensions of suitable sets of links (or sites) of the Ising 
percolation model at the critical point. In particular, y = 15/8 is the fractal dimension of the 
percolating cluster, y = 1 is related to the correlation length, ?/ = 13/24 is the fractal dimension 
of the red bonds (see |5^), y = 5/3 is the fractal dimension of the percolating cluster in the 



presence of a boundary (see ||60|), and y = 11/8 is the fractal dimension of the hull (see ||6l[] 
Unfortunately, the operator in which we are interested, being irrelevant, cannot be realized as 
a fractal dimension, but the coincidence of the other indices strongly supports the idea that it 
should also appear as subleading dimension of some suitably chosen set of links. 

Some theoretical justification of this remarkable coincidence of critical indices and fractal 
dimensions can be found in an interesting conjecture that was proposed for the first time in |^ 
and then discussed in detail in [59| and [3^. According to this conjecture, Ising percolation 
is described by the q ^ 1 limit of the tricritical g-state Potts model in exactly the same 
way in which the q 1 limit of ordinary g-state Potts describes standard percolation. The 
operator content of the q ^ 1 limit of the tricritical g-state Potts model can be studied with 
the same Coulomb gas techniques discussed above. It turns out that it contains (together 
with other operators) the nonunitary extension of the Ising model, and thus explains the above 
coincidence of critical indices and fractal dimensions. Notice that this conjecture is further 
supported by the identification as fractal dimensions of suitable sets of links of other critical 
indices that belong to the g — > 1 limit of the tricritical Potts but that are outside the nonunitary 
Ising class — see [BO] for a discussion. 



B Resummation of simple test functions 

In this appendix we consider a simple test function which behaves as (H) and whose perturbative 
expansion around g = is divergent but Borel summable. We show that many terms are needed 
in order to obtain the correct results, and, even worse, that in this case the standard method 
to set the error bars does not work properly. The estimated errors are much smaller than the 
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difference between the estimate and the exact value. 
Consider the function 



\ r+QO . , 

^^^^ " J-oc ~ ■ ^^-^^ 

Its expansion in powers of g, Z{g) = J2k ^ko'^, is Borel summable, and the large-order behavior 
of the fcth-order coefficient is given by 

Zk = oc {-l)\k - 1)! [1 + 0{l/k)] . (B.2) 

The function Z{g) is analytic in the complex plane with a cut along the negative real axis, and 
in particular it is analytic for (7 = 1 . For 6 = 1— g^Oit behaves as 

Z{g) = Zo + Zi6 + 0{6^), (B.3) 

where Zq = 0.9189189... and Zi = —0.0573155... In this case, in which the function is analytic, 
the resummation method we presented in the Introduction provides good estimates of the 
constants appearing in ( [B.3| ). One indeed obtains Zq = 0.9189(1) and Zi = —0.0572(3) from 
the 5th-order series, and Zq = 0.918919(1) and Zi = —0.057315(3) from the lOth-order series.^ 
Most important, the method provides correct estimates of the errors. 

In order to reproduce a nonanalytic behavior similar to ([21|), we consider the function 

B{b,g) = Z{g) + cil-gy+''. (B.4) 

Setting c = Zi, we have for (7 — 1 

5(6,(7) = Zo + Zi5 (1 + 5'') +0(5'). (B.5) 

We apply the same resummation procedure used for Z{g) to the perturbative expansion of 
B{b,g). To reproduce the correction predicted in the Ising case, we fix 6 = 1/7. The results of 
the analysis are now much less satisfactory. Indeed, we find Zq = 0.916(6) and Zi = —0.112(2) 
from the 5th-order series, and Zq = 0.918(1) and Zi = —0.103(6) from the lOth-order series. 
The estimate of Zq is not as precise as before, but the error is still correct. This is not surprising 
since the nonanalyticity is here rather weak, the nonanalytic corrections being of order 6^^^. 
On the other hand, the estimate of Zi, which is determined by resumming the dB{b, g)/dg 
(here, the nonanalytic corrections are stronger, of order 5^), is very imprecise and the estimate 
of the error, which is obtained from the stability analysis, is completely incorrect: the five- loop 
estimate differs from the exact value by more than 25 estimated error bars! Moreover, extending 
the series appears to be of little help. We conjecture that a similar phenomenon is happening 
in the FT estimates for = 1. Although the perturbative results indicate uJi ~ 4/3 with a 
tiny error, the correct result is sensibly different. 

^ The estimates and their errors are obtained using the procedure of . The estimate is obtained from the 
"optimal" values of the two free parameters introduced in the procedure (6 and a), which are determined by 
maximizing the stability of the results with respect to the order of the series analyzed. The errors are related to 
the stability of the results with respect to variations of the free parameters b and a around their optimal values. 
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We have also considered the case in which we add a term of the form Zig/ log(l — g), which 
mimicks the behaviour of the /3-function for > 3, observing completely analogous deviations. 

We have repeated the exercise by considering a nonanalytic singularity similar to that ex- 
pected in three dimensions, i.e. by setting 6 ~ 1. For example, for 6 = 9/10 we find Zq = 
0.917(1) and Zi = -0.068(4) from the 5th-order series, Zq = 0.9186(1) and Zi = -0.060(2) 
from the lOth-order series. As expected, the effect of the nonanalyticity is much smaller and 
the errors reasonable, although slightly underestimated. 



C Asymptotic expansion of large- integrals 

In this Appendix we wish to compute the asymptotic expansion for B ^ of integrals of the 
form 

where f{u) ~ for u — oo, and n and p are integers satisfying n > 1, p > 2. The function 
S{u) is given by 

6{u) = ~-\og\^ (C.2) 

where 



au) = J^ (C.3) 
V + 4 

The results presented here extend App. A of [|1^ to two dimensions. We wish to compute 
the leading nonanalytic contributions to the asymptotic expansion. For this purpose, we can 
replace 6{u) and f{u) with their leading behaviour for m ^ oo and write 

Ji/A [B + (2 log u) juy- 
where A is an arbitrary cutoff satisfying < A < 1. Then we make the substitution 

2 



u 



\ogu = y. (C.5) 



For y ^ 0, u oo, Eq. (|C.5D can be solved, obtaining the asymptotic expansion 

u- 21og(,/2)l'+„^,i/"- log"(,/2) /• ^^-^^ 

The first coefficients are: an = 022 = 1, ct2i = — 1- 

Substituting this expression in (|C.4|) and keeping only the leading contributions, we obtain 



/„(/,B)^ r t-wnl (^-^^ 
JO yuP ^ (B + y)" 



where analytic terms have been systematically neglected. 
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Since p > 2, we see that /«(/, ©) can be written as a sum of terms of the form 

Kr^mAO) - I dy + (C.8) 

with m, n, and p integers. The nonanalytic terms arc due to the integrals with n > 1, and thus 
we consider only this case. Now, observe that we need to consider n = 1 only , since 

Then, note also that 

logy)' 
e 

Thus, it is enough to consider Kiao, where a is a real number. In the following we assume 
a > 1. The final result however will be correct for all values of a. To compute the asymptotic 
expansion, first perform a Mellin transformation, rewriting 

-l/2+ioo ds 



[log(-logy)]^' = lim 



(C.IO) 



^^^e^ R^{A,s), (C.ll) 



where 



The previous equation defines Rai-A, s) for Re s < 0. By rotating the t contour one can obtain 
an analytic continuation in the domain Re s > with a cut along the positive real axis. In the 
following, we need the discontinuity at the cut. A simple calculation gives 

i?„(A, s+) - R^(A, S-) = f ^e^* = s^-\ (C.13) 

Jc t" 1 [a) 

where C is a contour running counterclockwise around the negative t-axis. We also need 
i?a(A, 0) = (— log A)^~"/(q; — 1). In order to compute the asymptotic expansion of Kiao{Q), 
deform the s-integral, so that it goes around the positive s-axis. Keeping into account the pole 
at s = we obtain 

c 77- r d ^ TT 

KiaoiQ) = Ra{^ 0) - / —-. — e^K(A, s+) - i?«(A, s-)] - / —-. — e^i?„(A, S), 

Jo 2m sm ITS Jc++C- 2msm'KS 

(C.14) 

where < < 1 is arbitrary and C± = {s : Re s — /i, ±Ims > 0}. The integral over the 
lines C± is of order and can therefore be discarded. In order to compute the integral over 
the cut, we make the substitution —s log Q = t, expand the integrand in powers of 1/ log and 
replace the upper integration limit — /ilog© with cxd — again we make an error of order G'^. 
The final integrations are trivial. We obtain finally 

i^i.o(e) ^ ^(-iogA)-« - (-ioge)^-"f:6 /^%+^~^^ ii^T^ (C-15) 

a-1 ^ r{a) \iogej 

where the coefficients bk are defined by 

= j:bkx''-\ (C.16) 

smx f-l 

fe=0 
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